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Abstract: This note generalizes the Fibonacci primitive roots to the set of integers. An asymp¬ 
totic formula for counting the number of integers with such primitive root is introduced here. 
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1 Introduction 

Let p > 5 be a prime. A primitive root r ^ ±1, modulo p is a Fibonacci primitive root if 
= r -|- (1 modp). Several basic properties of Fibonacci primitive roots were introduced in 
[38] . A few other plausible generalizations are available in the literature. Fibonacci primitive 
roots have applications in hnite held calculations, and other area of discrete mathematics, na, 
[H]. The complete statement appears in Theorem 8. 

The subset of primes with Fibonacci primitive roots is dehned by 

Vf = {p & ^ ordp(r) = (p(p) and r'^ = r + 1 mod p} C P. (1) 

The asymptotic formula for the counting function has the form 

Pf{x) = [p < X '. ordp(r) = (p(p) and = r -|- 1 mod p} (2) 

conditional on the generalized Riemann hypothesis, was established in [25] , and [39] . This note 
introduces a generalization to the subset of integers 

Mf = {n G M : ord„(r) = \{n) and = r + 1 mod n} C N, (3) 

and provides the corresponding asymptotic counting function 

Nf{x) = {n < X : ord„(r) = \{n) and = r -|- 1 mod n} . (4) 


Theorem 1. Assume the generalized Riemann hypothesis. Then, are infinitely many compos¬ 
ite integers n>l which have Fibonacci primitive roots. Moreover, the number of such integers 
n < X has the asymptotic formula 


Nf{x) 


filF—lOLF 

F (Oir) 





( 5 ) 


where A = {p G P : ordp(r) = p — 1, ordp 2 (r) ^ p{p — 1) or ^ r -|- 1 mod p^}, and Op = 
0.265705... is a constant. 
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The constant ap = (27/38) np> 2 (l ~ ^/p{p ~ 1)) is the average density of the densities of the 
primes in the residue classes mod 20 that have Fibonacci primitive roots, confer [38] for the 
calculations. 

Sections 2 to 7 provide some essential background and the required basic results. The proof of 
Theorem 1 is settled in Section 8. And Section 9 has a result for the harmonic sum constrained 
by Fibonacci primitive roots. 


2 Some Arithmetic Functions 

The Euler totient function counts the number of relatively prime integers (p(n) = ^{k : 
gcd(k,n) = 1}. This counting function is compactly expressed by the analytic formula <p(n) = 
^np|n(i - e N. 

Lemma 2. (Fermat-Euler) If a G Z is an integer such that gcd(a,n) = 1, then = 1 
mod n. 


The Carmichael function is basically a refinement of the Euler totient function to the finite ring 
Z/nZ. Given an integer n = p'fp'^f '' 'PTi Carmichael function is defined by 


A(n) = lcm(A(p);i),A(p2')---A(p/*)) = JJ (6) 

p”||A(n) 


where the symbol p“||n, z/ > 0, denotes the maximal prime power divisor of n > 1, and 


f (f (p'") if p > 3 or n < 2, 

I 2^“^ if p = 2 and v > 3. 


(7) 


The two functions coincide, that is, <p(n) = A(n) if n = 2,4,p™', or 2p™’,m > 1. And <p (2™) = 
2A (2”^). In a few other cases, there are some simple relationships between <p(n) and A(n). In 
fact, it seamlessly improves the Fermat-Euler Theorem: The improvement provides the least 
exponent A(n)|<p(?7,) such that = 1 mod n. 


Lemma 3. ([3 p. 233]) (i) For any given integer n G N the congruence = 1 mod n is 
satisfied by every integer a > 1 relatively prime to n, that is gcd(a,n) = 1. 

(ii) In every congruence = 1 mod n, a solution u exists which is a primitive root mod n, 
and for any such solution u, there are (p(A(n)) primitive roots congruent to powers of u. 


Proof: (i) The number A(n) is a multiple of every A (p^) such that p"\n. Ergo, for any relatively 
prime integer a > 2, the system of congruences 

= 1 mod pI\ = 1 mod p^^ ..., = 1 mod p/‘, (8) 


where t = uj{n) is the number of prime divisors in n, is valid. ■ 

An integer u G Z is called a primitive root mod n if min {m G M : u™' = 1 mod n} = A(n). 
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3 Characteristic Function For Fibonacci Primitive Roots 

The symbol ordpfc(M) denotes the order of an element r G Z)^ in the mnltiplicative 

group of the integers modulo p^. The order satisfies the divisibility condition ordpfe(r)|A(n), 
and primitive roots have maximal orders ordp*;(r) = A(n). The basic properties of primitive 
root are explicated in [2], [35], et cetera. The characteristic function / : N —> {0,1} of a hxed 
primitive root r in the hnite ring Z/p^ Z, the integers modulo p^, is determined here. 

Lemma 4. Letp^,k > 1, be a prime power, and let r G Z be an integer such that gcd (r,p^^ = 1. 
Then 


(i) The characteristic f function of the primitive root (r mod p^) is given by 


/(p‘) 



0 if p < 3, A; < 1, 

1 if ordpfc(r) = p^“^(p — 1) and = r + 1 mod p^, for any p > 5, k > 1, 

0 if ordpfc(r) ^ p^“^(p — 1), and p > 5, k > 2. 

(9) 


(ii) The function f is multiplicative, but not completely multiplicative since 
(hi) f{pq) = f{p)f{q), gcd(p, g) = 1, 

(iv) / (p^) ^ f {p)f {p), if ordp 2 (r) ^ p(p - 1). 

Proof: The congruences — x — 1 = 0 mod 2 and x^ — x — 1 = 0 mod 3 have no roots, so 
/ (2^) = 0 and / (3^) = 0 for all /c > 1. Moreover, for primes p > Sandfc > 1, function has the 
value / (p^) = 1 if and only if the element r G (Z /p^ Z) ^ is a Fibonacci primitive root modulo 
p^. Otherwise, it vanishes: / (p^) = 0. The completely multiplicative property fails for p = 5. 
Specihcally, 0 = / (5^) ^ /(5)/(5) = 1. ■ 

Observe that the conditions ordp(r) = p — 1 and ordp 2 (r) ^ p(p — 1) imply that the inte¬ 
ger r 7 ^ ±1, cannot be extended to a primitive root mod p^, k > 2. But that the condition 
ordp 2 (r) = p(p — 1) implies that the integer r can be extended to a primitive root mod p^,k > 2. 


4 Wirsing Formula 

This formula provides decompositions of some summatory multiplicative functions as products 
over the primes supports of the functions. This technique works well with certain multiplicative 
functions, which have supports on subsets of primes numbers of nonzero densities. 

Lemma 5. (jlTl p. 71]) Suppose that f : N — > C is a multiplicative function with the following 
properties. 

(i) f{n) > 0 for all integers n G N. 

(ii) / (p^) < for all integers k eN, and c < 2 constant. 

(hi) /(p) = {t + o(l))x/ logx, where t > b is a constant as x —> oo. 

p<x 

Then 
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n<x 


n] = 


eT'^r(r) 


+ 0 ( 1 ) 


X 


logo: 


n 

p<x 




P 


pZ 


( 10 ) 


The gamma function appearing in the above formula is defined by r(s) = s G C. 

The intricate proof of Wirsing formula appears in [47]. It is also assembled in various papers, 
such as mi, ESI P- 195], and discussed in RUl P- 70], ESI P- 3081. Various applications are 
provided in [27|, |S2], ES], et alii. 


5 Harmonic Sums And Products Over Fibonacci Primes 

The subset of primes Vp = {p G P : ord(r) = (p(p)andr^ = r + 1 mod p} C P consists of all 
the primes with Fibonacci primitive roots r G Z. By Theorem 8, it has nonzero density 
ap = S iVp) > 0. The real number Oi? > 0 is a rational multiple of the well known Artin 
constant, see [SS], and for the derivation. 

The proof of the next result is based on standard analytic number theory methods in the liter¬ 
ature, refer to [SSI Lemma 4]. 

Lemma 6. Assume the generalized Riemann hypothesis, and let x > 1 be a large number. 
Then, there exists a pair of constants (3p > 0, and > 0 such that 


(i) 


E - 

p<x,p&Vp 


ap log logx + /3p + O 


/log logx\ 

V log^ X J ■ 


(ii) 


E 

P<x,p£Vf 


logp 

p — 1 


ap logx — 'yp + O 


/ log log X \ 
V ^og^x )' 


Proof: (i) Let 'xp{x) = {p < x : ord(r) = (p(p) and = r -|- 1 mod p} be the counting mea¬ 
sure of the corresponding subset of primes Vp. To estimate the asymptotic order of the prime 
harmonic sum, use the Stieltjes integral representation: 




TTp{x) 


p<x,pe'PF 


X 


+ C(Xo) + 


fXQ 




dt, 


( 11 ) 


where Xq > 0 is a constant. By the GRH, Tip{x) = ap7i{x) = apx/ logx + O (log logx/log^ x), 
see Theorem 8. This yields 



ap ^ ^ / log logx \ 
logx V log^(x) J 


+ Co(a^o) 


Fap 


I XQ 


+ 0 


/ log logt 
\tlog\t) 


dt 


Jlogt 

ap log log X - log log xo + Co (xq) + O 


log log x\ 

log^ X / 


( 12 ) 


where (3p = — log log Xq + Cq (xq) is a generalized Mertens constant. The statement (ii) follows 
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from statement (i) and partial snmmation. ■ 

A generalized Mertens constant and a generalized Enler constant 7 ^? have other eqnivalent 
dehnitions snch as 


^ log logx j and = 7u - (13) 

\p<x,pGVp ^ / P&'Pf, k>2 ^ 

respectively. These constants satisfy Pf = (diap and 7 ^? = 707 ?. If the density ap = 1, these 
dehnitions rednce to the nsnal Euler constant and the Mertens constant, which are dehned by 
the limits 


7 = 




and j3i 



(14) 


or some other equivalent dehnitions, respectively. Moreover, the linear independence relation 
in flT^ becomes /9 = 7 — X)p >2 Z)fc >2 {kp’^) see [El Theorem 427]. 


A numerical experiment utilizing a small subset of primes with Fibonacci primitive roots gives 
the followings approximate values: 




p>2 


P{P - 1 ) 


= 0.265705484288843681890137..., 


ii) V --ap log log a; = 0.05020530308647012230491, and 

77 


p<1301,p€7^F 


p 


Uii) 7 f 


E 


p<130l,p£'Pp 


logp 

P — 1 


- ap logx = 0.0221594862523476326826286 .. 


Lemma 7. Assume the generalized Riemann hypothesis, and let x > 1 be a large number. 
Then, there exists a pair of eonstants 7 ^ > 0 and up > 0 sueh that 


n b 


p<x,p£'Pp 


1 \ ^ _ , _ „ /log log X 


p 


= log(x)"^ + O 


V log^ 


X 


p<x,pCiPp 


P 


(ii) n + Yl i^~P log(a^)"^ + O 


P&Vp 


-1 


/ log log X 
V log^ X 

X°^F iQg iQg ^ 


(iii) TT (l - = e*''’-’"!”' +0 . , 

vAhvA V log X 

Proof: (i) Express the logarithm of the product as 

-1 


E '°s(i-p 

P<x,pG'Pf ^ 


E E/f" E E Ei 

P<x,pG'Pf, k>l p<x,p^'Pf P^x^pGVfi k>2 


(16) 


Apply Lemma 6 and relation ffT3D to complete the verihcation. For (ii) and (iii), use similar 
methods as in the hrst one. ■ 
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The constant > 0 is defined by the donble power series (an approximate nnmerical valne for 
p < X = 1301 is shown): 


1/2 = 



^ 0.188622600886988493134287... . 


( 16 ) 


6 Density Correction Factor 

The snbsets of primes Vp = {p G P : ord(r) = p — 1 and = r + 1} has an important partition 
as a disjoint nnion Vp = AUB, where 

(i) A= [p : ord(r) = p — 1, ordpkV ^ — 1), or ^ r + 1 mod and 

(ii) B = [p gF : ord(r) = p — 1, ordp^r = p^“^(p — 1), and = r + 1 mod p^} . 

This partition has a role in the determination of the density of the snbset of integers Afp = 
{n G M : ord(r) = X{n) and = r + 1 mod n} with Fibonacci primitive roots r G Z. The finite 
prime prodnct occnrring in the proof of Theorem 1 splits into two finite prodncts over these snb¬ 
sets of primes. These are nsed to reformnlate an eqnivalent expression snitable for this analysis: 


P{x) 


n 

ordp(r)=p—1, 

ordp 2 1) or r^^r+lmod 

n 

p^<x, 

ordp(r)=p-l, 

oi’dp 2 (r) 7 ^p(p—1) or r^^r+lmod 


n b-4 n b 


1 

1 H— 

p 


1 - 


pz 


P<x,pGA 


P 


P<x,p&'Pf 


p£A 


pz 


nb-b n b4 +« 


P^XyP^Vp 


P 


n 

p^<x, 

ordp 2 (t)=p(p— 1) and r^=r+lmod p^ 


n 

p^<x, 

ordp 2 (r)=p(p—1) and r^=r-|-lmod p^ 


-1 


1 1 

1 H-1—5- -l- ■ ■ ■ 

p p^ 


1 - 


p 


-1 


(17) 


logx 


X 


where the convergent partial prodnct is replaced with the approximation 



( 18 ) 


The prodnct Hpe/t ~ P~‘^) rednces the density to compensate for those primes for which 
the Fibonacii primitive root r mod p that cannot be extended to a Fibonacci primitive root 
r modp^. This seems to be a density correction factor similar to case for primitive roots over 
the prime nnmbers. The correction reqnired for certain densities of primes with respect to fixed 
primitive roots over the primes was discovered by the Lehmers, see BB. 
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7 The Proof Of The Theorem 

The algebraic constraint = r + 1 restricts the Fibonacci primitive roots r modnlo a prime 
p > 5 to the algebraic integers (l ± \/5) / 2 in the finite field Fp. This in tnrns forces the 
corresponding primes to the residne classes p = ±1 mod 10, the analysis appears in [38]. A 
conditional proof for the conjectnre of Shank was achieved in [25| and [39]. 

Theorem 8. ([25]) // extended Riemann hypothesis hold for the Dedekind zeta function 
over Galois fields of the type Q v^), where n is a squarefree integer, and d\n. Then, the 
density of primes which have Fibonacci primitive roots is ap = (27/38) np> 2 (l “ l/p(p “ 1))^ 
and the number of such primes is 


P(r)-n ^ ^fxloglogx 

1 pyX j — 0:p -“T cy I - 9 - I as X 

log X \ log X 


oo. 


(19) 


Proof of Theorem 1: By Theorem 8, the density ap = 6 ifPp) > 0 of the snbset of primes 
Vp is nonzero. Pnt t = ap in Wirsing formnla, Lemma 6, and replace the characteristic fnnc- 
tion f{n) of Fibonacci primitive roots in the finite ring Z Z, /c > 1, see Lemma 4, to prodnce 


n<x 


n] = 


1 X 

+ o(l) 


eT'^F(r 

1 

eT'“^'F {ap) 


logx 


+ 0 ( 1 ) 


p<x 

X 




logo: 


f{p) f{p^) 

9 

P 

n 


+ 


p’^<x, 

ordp(r’)=p— 1 , 

oi'd„2 (i’)7^p(p— 1 ) or r^^r+lmod 


1 + 


p 


X 


n 

ord 2('P)=p(p~l) r^=r+lmod 


1 1 

1 H- 1 — 5 - + 

p p^ 


( 20 ) 


Note: In terms of the snbsets A and B, see Section 6, the characteristic fnnction f{n) of 
Fibonacci primitive roots has the simpler description 

f / 1 if p e and fc = 1, , ^ ^ r 1 if peB, 

J \P ) i g if p G 9l, and k >2, 7 [P ) S g otherwise. 


Replacing the eqnivalent prodnct, see ffT7|l in Section 6, and nsing Lemma 7, yield 




n<x 


e'y^FV {ap) 


+ 0 ( 1 ) 


X 


logx 


peA 


nb-^ n 




F {ap) 


+ 0 ( 1 ) 


X 


(logo:) 


i-ap n 


p€A 


pz 


pz 


P<x,p£'Pf 


1 - 


p 


-1 


( 22 ) 


where Oi? > 0 is a constant, and jp is a generalized Enler constant, see Lemma 6 for detail. 


The constant ap = (27/38) np> 2 (l “ 7/p{p — 1)) = 0.265705... is the average density of the 
densities of the primes in the residne classes mod 20 that have Fibonacci primitive roots, confer 
[38] for the calcnlations. 
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8 Harmonic Sum For Fibonacci Primitive Roots 

The first few primes in the subset of primes Vf = {p < x : ord(r) = (p(p) and = r + 1 mod p}, 
which have Fibonacci primitive roots, are 

Vf = {5,11,19, 31,41, 59, 61, 71, 79,109,131,149,179,191, 239, 241, 251, 269, 271, 311, 

359, 379, 389,409,419,431,439,449,479,491,499, 569, 571, 599, 601,631,641,659, 
701, 719, 739, 751, 821,839, 929, 971,1019,1039,1051,1091,1129,1171,1181, 
1201,1259,1301,...}. (23) 

The subset of primes Vf C P of density ap = S(Vf) > 0 generates a subset of composite 
integers A/f = {n G N : ord(r) = A(?7,) and = r + 1 mod n}, which have Fibonacci primitive 
roots. The hrst few are listed as 

^^F = {5,11,19, 31,41, 55, 59, 61, 71, 79, 5 • 19,109, ll^, 5 • 19,149, 5 • 31,...}. (24) 


The hrst such composite number n = 55 has two Fibonacci primitive roots r = 8 and 48 of 
order ordr = A(55) = 20. Let 


Mf = {n G N : ord(r) = A(n) and = r + 1 mod n} , 


and let 


Nf{x) = [n < X : ord(r) = A(n) and = r + 1 mod n} 
be the corresponding discrete counting measure. This was determined in Theorem 1. 


(25) 

(26) 


The subset Mf C 77 is a proper subset TZ = {n E N : p\n p E Vf}- This sub¬ 
set of integers, which contains all the primes powers, is generated by the subset of primes 
= {p G P : ord(r) = A(p) and r‘^ = r + 1 mod p}. Since there are some primitive root 
(r mod p) which are not primitive root modulo prime powers p™, m >2, the subset 77 is slightly 
larger than the subset F/f- For example, r = 3 is a Fibonacci primitive root mod p = 5, but 
(3 -|- 5fc)^ ^ 3 + 5k mod 5^ for k E Fp. Thus it cannot be extended to a Fibonacci primitive 
root modulo 5™, m > 2. More precisely, the prime powers 5™ G 77, m > 1, but 5™ ^ AfF, m >2. 

An asymptotic formula for the harmonic sum over the subset of integers TVf is determined here. 


Lemma 9. Let x > 1 be a large number, let aF = S {Vf) > 0 be the density of the subset of 
primes Vf, and let Mf gN be a subset of integers generated by Vf- Then 


^ - = KF(logx)"^+7F + 0 


rtMT nG- Kf 1 


nr^cr l—otF 


(27) 


The number of = 0.265705484288...is a generalized Artin constant, and > 0 is a generalized 
Euler constant, see 0131) for the dehnition. The other constant is 


Kf 


a ^F (cK f) 



( 28 ) 


where the index of the product ranges over the subset primes A, see Section 6. 
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Proof: Use the discrete counting measure Nf{x) = [apSiF + o(l)) x(logx)"^ Theorem 1, to 
write the hnite sum as an integral, and evaluate it: 


E - 

^ . n 

n<x,n&NF 



Npit) 

t 


X 


+ 


XQ 



Npit) 

F 


dt, 


(29) 


where Xq > 0 is a constant. Continuing the evaluation yields 



(cufK-f + c»(l)) 


+ Co (xo) + 


(logx)^“"^ 
KFilogx)'^^ +-fF + 0 


{apf^F + 0 ( 1 )) 

Lo t{log 

1 


dt 


(logx)^”"-?' 


(30) 


where Cq (xq) is a constant. Moreover, 


■jF = lim 

x—^oo 



Co (xo) + 


/’°° {api^F + 

Jxo yiogty-^^ 


is a second definition of this constant. ■ 


(31) 


The integral lower limit xq = 4 appears to be correct one since the subset of integers is 
AA 2 = {5,11,...}. 

The Fibonacci primitive root problem is a special case of a more general problem that inves¬ 
tigates the existence of infinitely many primes such that the polynomial congruence f{n) = 
0 mod p has primitive root solutions. The algebraic and analytic analysis for these more general 
polynomials are expected to be interesting. 


9 Problems 

Problem 8.1 Let A = {p ElP : ord(r) = p — 1, ord(r) y p{p — 1) or ^ r -|- 1 mod p^}. De¬ 
termine whether or not the subset of primes M = { 5, ... } is finite. There is an algebraic 
characterization, see (16). Sketch of the proof: Use a) = r -|- 1 mod p, b) 1 < r < p to show 
that (r -|- mpY = r + mp + mod p^, m G Fp has no solutions for all primes p > 5. 

Problem 8.2 Let Lp{n)/\{n) be the index of the finite group 'L/n'L, and suppose it contains 
a Fibonacci primitive root. Explain which maximal subgroup G C. (Z/p^ of cardinality 
y^G = \{n) contains the Fibonacci primitive root. Some theoretical/numerical data is available 

in [7j. 

Problem 8.3 Show that the constants satisfy the equation 7^7 — 70 ^? = 0, confer equation ([13]). 
Problem 8.4 Determine the least primitive fibonacci root modulo n. 

Problem 8.5 Show that L{s,xf) = 'n,n>iXF{n)/n^ is analytic on the half palne 3fJe(s) > 1 , 
and has a pole at s = 1. Here the character is defined by = 1 if n is a fibonacci integer, 

else it vanishes. 
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Problem 8.6 Show that the logarithm derivative —L'(s,xf)/L(s,xf)^ = X]n>i 

analytic on the half palne 3fJe(s) > 1, and has a pole at s = 1. Here the prime characteristic 

fnnction is dehned by Aj^(n) = logn if n is a hbonacci prime, else it vanishes. 

Problem 8.7 Generalize the analysis for Fibonacci primitive root solntions of the qnadratic 
eqnation — x — 1 = 0 mod n for inhnitely many integers n > 1 to the existence of primitive 
root solntions for a general quadratic equation ax"^ + bx + c = 0 mod n for inhnitely many 
integers n > 1. 
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